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Analysis of the radio-metric tracking data from the Pioneer 10 and 11 spacecraft at distances 
between 20 and 70 astronomical units from the Sun has consistently indicated the presence of 
an anomalous, small, and constant Doppler frequency drift. The drift is a blueshift, uniformly 
changing at the rate of (5.99 ± 0.01) x 10 -9 Hz/s. The signal also can be interpreted as a constant 
acceleration of each particular spacecraft of (8.74 ± 1.33) x 10 _8 cm/s 2 directed toward the Sun. 
This interpretation has become known as the Pioneer anomaly. We provide a problem set based on 
the detailed investigation of this anomaly, the nature of which remains unexplained. 

PACS numbers: 04.80.-y, 95.10.Eg, 95.55.Pe 



I. INTRODUCTION 

The Pioneer 10/11 missions, launched on 2 March 1972 
(Pioneer 10) and 4 December 1973 (Pioneer 11), were the 
first to explore the outer solar system. After Jupiter and 
(for Pioneer 11) Saturn encounters, the two spacecraft 
followed escape hyperbolic orbits near the plane of the 
ecliptic to opposite sides of the solar system. (The hy- 
perbolic escape velocities are ^12.25 km/s for Pioneer 10 
and ~11.6 km/s for Pioneer 11.) Pioneer 10 eventually 
became the first man-made object to leave the solar sys- 
tem. See Fig. 1 for a perspective of the orbits of the 
spacecraft and Tables I and II for more details on their 
missions. The orbital parameters for the crafts are given 
in Table III. 

The Pioneers were excellent craft with which to per- 
form precise celestial mechanics experiments due to a 
combination of many factors, including their attitude 
control (spin-stabilized, with a minimum number of at- 
titude correction maneuvers using thrusters), power de- 
sign [the Plutonium-238 powered heat-source radioiso- 
tope thermoelectric generators (RTGs) are on extended 
booms to aid the stability of the craft and also reduce 
the effects due to heating], and precise Doppler tracking 
(with the accuracy of post-fit Doppler residuals at the 
level of mHz). The result was the most precise naviga- 
tion in deep space to date. See Fig. 2 for a design drawing 
of the spacecraft. 

By 1980 Pioneer 10 had passed a distance of ~20 as- 
tronomical units (AU) from the Sun and the accelera- 
tion contribution from solar-radiation pressure on the 
craft (directed away from the Sun) had decreased to 
< 4 x 10 _8 cm/s 2 . At that time an anomaly in the 
Doppler signal became evident and persisted. 1 Subse- 
quent analysis of the radio-metric tracking data from the 
Pioneer 10/11 spacecraft at distances between ~20 and 
70 AU from the Sun has consistently indicated the pres- 
ence of an anomalous, small, and constant Doppler fre- 



quency drift of v = (5.99 ± 0.01) x 10" 9 Hz/s. The drift 
can be interpreted as being due to a constant accelera- 
tion of the spacecraft of a P = (8.74 ± 1.33) x 10~ 8 cm/s 2 
directed toward the Sun. 2 The nature of this anomaly 
remains unknown, and no satisfactory explanation of the 
anomalous signal has been found. This signal has become 
known as the Pioneer anomaly. 3 

Although the most obvious explanation would be that 
there is a systematic origin to the effect, perhaps gen- 
erated by the spacecraft themselves from excessive heat 
or propulsion gas leaks, 4-6 none has been found. 7-9 Our 
inability to explain the anomalous acceleration of the Pi- 
oneer spacecraft has contributed to a growing discussion 
about its origin, especially because an independent veri- 
fication of the anomaly's existence has been performed. 10 

Attempts to verify the anomaly using other spacecraft 
have proven disappointing, 2:9:11,12 because the Voyager, 
Galileo, Ulysses, and Cassini spacecraft navigation data 
all have their own individual difficulties for use as an in- 
dependent test of the anomaly. In addition, many of the 
deep space missions that are currently being planned ei- 
ther will not provide the needed navigational accuracy 
and trajectory stability of under 10~ 8 cm/s 2 (that is, 
recent proposals for the Pluto Express and Interstellar 
Probe missions), or else they will have significant on- 
board effects that mask the anomaly (the Jupiter Icy 
Moons Orbiter mission). The acceleration regime in 
which the anomaly was observed diminishes the value 
of using modern disturbance compensation systems for 
a test. For example, the systems that are currently be- 
ing developed for the LISA (Laser Interferometric Space 
Antenna) and LISA Pathfinder missions are designed to 
operate in the presence of a very low frequency accel- 
eration noise (at the millihcrtz level), while the Pioneer 
anomalous acceleration is a strong constant bias in the 
Doppler frequency data. In addition, currently available 
accelerometers are a few orders of magnitude less sensi- 
tive than is needed for a test. 11-15 Should the anomaly 
be a fictitious force that universally affects frequency 
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FIG. 1: Ecliptic pole view of the Pioneer 10, Pioneer 11, and Voyager trajectories. Pioneer 10 is traveling in a direction almost 
opposite to the galactic center, while Pioneer 11 is heading approximately in the shortest direction to the heliopause. The 
direction of the solar system's motion in the galaxy is approximately towards the top. (Digital artwork by T. Esposito, NASA 
ARC Image # AC97-0036-3.) 



TABLE I: Pioneer 10/11 mission status. 





Pioneer 10 


Pioneer 11 


Launch 


2 March 1972 


5 April 1973 


Planetary encounters 


Jupiter: 4 Dec 1973 


Jupiter: 2 Dec 1974 
Saturn: 1 Sep 1979 


Last data received 
at heliocentric distance 


27 April 2002 
w 80.2 AU 


1 October 1990 
« 30 AU 


Direction of motion 
time to reach destination 


Star Aldebaran 
w 2 million years 


Constell. of Aquila 
~ 4 million years 



standards, 2 the use of accelerometers will shed no light 
on what is the true nature of the observed anomaly. 

A comprehensive test of the anomaly requires an es- 
cape hyperbolic trajectory, 2 ' 9 ' 11 which makes a number 
of advanced missions (such as LISA, the Space Test of 
Equivalence Principle, and LISA Pathfinder) less able to 
test properties of the detected anomalous acceleration. 
Although these missions all have excellent scientific goals 



and technologies, they will be in a less advantageous po- 
sition to conduct a precise test of the detected anomaly 
because of their orbits. These analyses of the capabili- 
ties of other spacecraft currently in operation or in design 
lead to the conclusion that none could perform an inde- 
pendent verification of the anomaly and that a dedicated 
test is needed. 11-14 Because the Pioneer anomaly has re- 
mained unexplained, interest in it has grown. 13-15 Many 
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TABLE II: Position of Pioneer 10 on 1 January 2005 (data 
taken from the Jet Propulsion Laboratory's Horizons system 
<http : //ssd . jpl . nasa . gov/horizons . html>) . 



Distance from Sun 


87.06 AU 


Position (lat., long.) 


(3.0°, 78.0°) 


Speed relative to the Sun 


12.24 km/s 


Distance from Earth 


13.14 x 10' J km 


Round- Trip Light Time 


« 24 hr 21 min 



TABLE III: Orbital parameters for Pioneers 10 and 11 at 
epoch 1 January 1987, 01:00:00 UTC. The semi-major axis is 
a, e is the eccentricity, / is the inclination, O is the longitude of 
the ascending node, u> is the argument of the perihelion, Mo is 
the mean anomaly, /o is the true anomaly at epoch, and ro is 
the heliocentric radius at the epoch. The direction cosines of 
the spacecraft position for the axes used are (a, (3, 7). These 
direction cosines and angles are with respect to the mean 
equator and equinox of J2000. The ecliptic longitude £0 and 
latitude 60 are also listed for an obliquity of 23°26 21. 4119. 
The numbers in parentheses denote realistic standard errors 
in the last digits. 



Parameter 


Pioneer 10 


Pioneer 11 


a [km] 


-1033394633(4) 


-1218489295(133) 


e 


1.733593601(88) 


2.147933251(282) 


/ [deg] 


26.2488696(24) 


9.4685573(140) 


fl [deg] 


-3.3757430(256) 


35.5703012(799) 


uj [deg] 


-38.1163776(231) 


-221.2840619(773) 


M [deg] 


259.2519477(12) 


109.8717438(231) 


fo [deg] 


112.1548376(3) 


81.5877236(50) 


r [km] 


5985144906(22) 


3350363070(598) 


a 


0.3252905546(4) 


-0.2491819783(41) 


13 


0.8446147582(66) 


-0.9625930916(22) 


7 


0.4252199023(133) 


-0.1064090300(223) 


£ [deg] 


70.98784378(2) 


-105.06917250(31) 


b [deg] 


3.10485024(85) 


16.57492890(127) 



novel ideas have been proposed to explain the anomaly, 
such as a modification of general relativity, 16-18 a change 
in the properties of light, 19 or a drag from dark matter. 20 
So far, none has been able to unambiguously explain the 
anomaly. 

The aim of the following problems is to excite and en- 
gage students and to show them by example how research 
really works. As a byproduct wc hope to demonstrate 
that they can take part and understand. 

In Sec. II we begin with a few calculations that are 
designed to provide intuition into the dynamics of the 
problem. Then we concentrate on the error analysis that 
demonstrated the existence of the anomaly in the data. 
In Sec. Ill we will consider possible sources of errors in 
the measurement of the Pioneer anomalous acceleration 
that have origins external to the spacecraft. Section IV 
will consider likely sources of errors that originate on- 
board of the spacecraft. We will discuss the anticipated 



effect of computational errors in Sec. V and summarize 
the various contributions in Sec. VI. In Appendix A we 
will discuss the Allan deviation - a quantity useful in 
evaluating the stability of the reference frequency stan- 
dards used for spacecraft navigation. 

II. EFFECT OF THE ANOMALY ON THE 
PIONEERS' TRAJECTORIES 

For these problems we can use the approximation of 
ignoring the angular momentum in the hyperbolic orbits, 
and treat the velocities as radial. 

Problem 2.1. Given the values of the orbital param- 
eters of Pioneer 10 on 1 January 2005 (distance r and 
velocity vo) in Table II, what would be the final velocity 
of Pioneer 10 assuming that there is no anomaly? 

Solution 2.1. The terminal escape velocity can be 
calculated from conservation of energy: 

Voo = ^ 2 -2GM Q /r = 11.38 km/s. (1) 

Problem 2.2. Assume that the anomaly continues 
and is a constant. At what distance from the Sun would 
the acceleration of the anomaly be equal to that of grav- 
ity? 

Solution 2.2. We equate the anomalous and Newto- 
nian gravitational forces and find 

d = \J GMq I a p w 2, 600 AU. (2) 

Problem 2.3. Assume the anomaly is due to a phys- 
ical force. The time interval over which the Pioneer 10 
data was analyzed was 11.5years (1987 to 1998.5). Be- 
cause of the anomaly, what is the approximate shortfall 
in distance which Pioneer 10 traveled during this time? 

Solution 2.3. Simple mechanics yields 

Ad = -apt 2 w 57, 500 km. (3) 

Problem 2.4. If the anomaly were to continue as a 
constant out to deep space and be due to a force, how far 
would Pioneer 10 travel before it reaches zero velocity and 
starts to fall back in to the Sun? When will the Pioneer 
reach zero velocity? 

Solution 2.4. As seen in Problem 2.1, the gravi- 
tational potential energy in deep space due to the Sun 
is small compared to the kinetic energy of the Pioneer. 
Therefore, we can ignore it compared to the kinetic en- 
ergy needed to be turned into potential energy. (This 
approximation slightly overestimates the time and dis- 
tance.) If we use the data from Table II, we can find the 
time by solving 

Vf = = vo — apt. (4) 

The solution is t = 1.40 x 10 13 s = 445, 000 yr. The dis- 
tance traveled would be given by 

D = r Q + v t-\a P t 2 w 573, 300 AU = 9 light years. (5) 
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FIG. 2: Drawing of the Pioneer spacecraft. 



This distance would be well on the Pioneer 10's way to 
Aldebaran (see Table I). Because this distance is very 
large compared to r = 87.06 AU, this result verifies 
the validity of the approximation of ignoring the grav- 
itational potential energy. 

Problem 2.5. What if the Pioneer anomaly is caused 
by heat from the RTGs? As stated, they are powered 
by Plutonium 238, which has a half life of 87.74 years. 
(See Sec. Ill for more details on the RTGs.) What would 
be the change in the final escape velocity of the Pioneers 
from this cause? 

Solution 2.5. The acceleration would be decreasing 
exponentially due to the radioactive decay of the 238 Pu. 
Therefore, the change in the velocity would be 

/>oo 

Av = a P / dt exp[-i(ln 2)/(87.74 yr)] 
Jo 

= a P (87.74 yr)/ln 2 = 3.5 m/s. (6) 

III. SOURCES OF ACCELERATION NOISE 
EXTERNAL TO THE CRAFT 

External forces can contribute to all three vector com- 
ponents of spacecraft acceleration. However, for a rotat- 
ing spacecraft, such as the Pioneers, the forces generated 



on board produce an averaged contribution only along 
the direction of the axis of rotation. The difference in 
the second case is due to the fact that the two non-radial 
components (those that are effectively perpendicular to 
the spacecraft spin) are averaged out by the rotation of 
the spacecraft. Furthermore, non-radial spacecraft accel- 
erations are difficult to observe by the Dopplcr technique, 
which measures the velocity along the Earth-spacecraft 
line of sight. Although we could in principle set up com- 
plicated models to predict effects due to all or each of the 
sources of navigational error, often the uncertainty of the 
models is too large to make them useful, despite the sig- 
nificant effort required to use them. A different approach 
is to accept our ignorance about these non-gravitational 
accelerations and assess to what extent they could con- 
tribute to spacecraft acceleration over the time scales of 
all or part of the missions. 

In this section we will discuss possible sources of accel- 
eration noise originated externally to the spacecraft that 
might significantly affect the navigation of the Pioneers 
10 and 11. We first consider the forces that affect the 
spacecraft motion, such as those due to solar-radiation 
pressure and solar wind pressure. We then discuss the 
effects on the propagation of the radio signal that are 
from the solar corona, electromagnetic forces, and the 
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phase stability of the reference atomic clocks. 

Problem 3.1. There is an exchange of momentum 
when solar photons impact the spacecraft and are either 
absorbed or reflected. The models of this process take 
into account various parts of the spacecraft surface ex- 
posed to solar radiation, primarily the high-gain antenna, 
predict an acceleration directed away from the Sun as a 
function of the orientation of the spacecraft and its dis- 
tance to the Sun. The effect of the solar radiation pres- 
sure depends on the optical properties of the spacecraft 
surface (that is, the absorptivities, reflectivities and emis- 
sivities of various materials used to build the spacecraft), 
and the effective areas of the spacecraft part. The effect 
can be distinguished from the 1/r 2 law due to gravity 
because the direction between the Sun and the effective 
spacecraft surface varies. 

For navigational purposes, we determine the magni- 
tude of the solar-pressure acceleration at various orien- 
tations using Doppler data. The following equation pro- 
vides a good model for analyzing the effect of solar radi- 
ation pressure on the motion of distant spacecraft, ft is 
included by most of the programs used for orbit determi- 
nation: 



a sp (r) 



nfo A cos 9 (r) 



(7) 



where /© = 1367 W/m 2 (AU) 2 is the effective Stefan- 
Boltzmann temperature or solar radiation constant at 
1 AU from the Sun, A is the effective area of the craft as 
seen by the Sun, and c is the speed of light. The angle 
6 is the angle between the axis of the antenna and the 
direction of the Sun, m is the mass of the spacecraft, and 
r is the distance from the Sun to the spacecraft in AU. 

For the Pioneers the effective area of the spacecraft, 
A, can be taken to be the antenna dish of radius 1.37 m. 
(See Fig. 2 for more information.) The quantity m is the 
mass of the Pioneers when half of the hydrazine thruster 
fuel is gone (241kg). Finally, n in Eq. (7) is the effective 
absorption/reflection coefficient of the spacecraft surface, 
which for Pioneer 10 was measured to be k = 1.71. A 
similar value was obtained for Pioneer 11. 

Estimate the systematic error from solar-radiation 
pressure on the Pioneer 10 spacecraft over the interval 
from 40 to 70.5 AU, and for Pioneer 11 from 22.4 to 
31.7 AU. 

Solution 3.1. With the help of Eq. (7) we can esti- 
mate that when the craft reached 10 AU, the solar radi- 
ation acceleration was 18.9 x 10~ 8 cm/s 2 decreasing to 
0.39 x 10 _8 cm/s 2 for 70 AU. Because this contribution 
falls off as r~ 2 , it can bias the Doppler determination 
of a constant acceleration. By taking the average of the 
r~ 2 acceleration curves over the Pioneer distance, we can 
estimate the error in the acceleration of the spacecraft 
due to solar-radiation pressure. This error, in units of 
10" 8 cm/s 2 , is 0.001 for Pioneer 10 over the interval from 
40 to 70.5 AU, and 0.006 for Pioneer 11 over the interval 
from 22.4 to 31.7AU. 

Problem 3.2. Estimate the effect of the solar wind 



on the Pioneer spacecraft. How significant is it to the 
accuracy of the Pioneers' navigation? 

Solution 3.2. The acceleration caused by the solar 
wind has the same form as Eq. (7), with /©/c replaced 
by m p v 2 n, where m p is the mass of proton, n k, 5cm~ 3 
is the proton density at 1 AU, and v 400km/s is the 
speed of the wind. Thus, 



r(r) 



nnipV 2 n Acos9(r) 
m r 2 

4.4 x 10" 



11 /20AUx2 2 
{—^— ) cm /s 



(8) 



We will use the notation a to denote accelerations that 
contribute only to the error in determining the magnitude 
of the anomaly. 

Because the density of the solar wind can change by 
as much as 100%, the exact acceleration is unpredictable. 
Even if we make the conservative assumption that the so- 
lar wind contributes only 100 times less force than solar 
radiation, its smaller contribution is completely negligi- 
ble. 

Problem 3.3. Radio observations in the solar system 
arc affected by the electron density in the outer solar 
corona. Both the electron density and density gradient 
in the solar atmosphere influence the propagation of radio 
waves through this medium. Thus, the Pioneers' Doppler 
radio-observations were affected by the electron density 
in the interplanetary medium and outer solar corona. In 
particular, the one way time delay associated with a plane 
wave passing through the solar corona is obtained by in- 
tegrating the group velocity along the ray's path, £: 



At = -- 



2cn crit 



1 f sc 
— TT / ' 

;ritW J© 



d£n e (t,r), 



(9) 



n clit (u) - 1.240 x lO^^^-) 2 cm" 3 , (10) 

where n e (i,r) is the free electron density in the solar 
plasma, and n cr i t {v) is the critical plasma density for the 
radio carrier frequency v. 

We see that to account for the plasma contribution, 
we should know the electron density along the path. We 
usually write the electron density, n e , as a static, steady- 
state part, n e (r), plus a fluctuation 5n e (t 7 r): 



n e (t, r) = n e (r) + 5n e (t, r). 



(11) 



The physical properties of the second term are diffi- 
cult to quantify. But its effect on Doppler observables 
and, therefore, on our results is small. On the con- 
trary, the steady-state corona behavior is reasonably well 
known and several plasma models can be found in the 
literature. 2 ' 21 ^ 24 

To study the effect of a systematic error from propaga- 
tion of the Pioneer carrier wave through the solar plasma, 
we adopted the following model for the electron density 
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profile, 



n e (t,r) = 



r 



where r is the heliocentric distance to the immediate ray 
trajectory and <j) is the helio-latitude normalized by the 
reference latitude of <po = 10°. The parameters r and 
</> are determined from the signal propagation trajectory. 
The parameters A, B, C describe the solar electron den- 
sity. (They are commonly given in cm -3 .) 

Define Ac? = cAt to be the delay of radio signals due 
to the solar corona contribution and use Eqs. (9)-(12) to 
derive an analytical expression for this quantity. 

Solution 3.3. The substitution of Eq. (12) into 
Eq. (9) results in the following steady-state solar corona 
contribution to the earth-spacecraft range model that 
was used in the Pioneer analysis: 2 



Ad 



+ 



D 



(13) 



Here r/ = Rq/2 n C rit(^o) is the units conversion factor 
from cm~ 3 to meters, vq is the reference frequency (y$ = 
2295 MHz for Pioneer 10), v is the actual frequency of the 
radio wave, and p is the impact parameter with respect 
to the Sun. The function F in Eq. (13) is a light-time 
correction factor that was obtained during integration of 
Eqs. (9) and (12) as: 2 



F = F(p, tt, Te) — —{ arctan 

7T 



arctan 



(14) 



where tt and te are the heliocentric radial distances to 
the target spacecraft and to the Earth, respectively. 

Problem 3.4. The analyses of the Pioneer anomaly 2 
took the values for A, B, and C obtained from the Cassini 
mission on its way to Saturn and used them as inputs 
for orbit determination purposes: (A, B, C)r] = (6.0 x 
10 3 , 2.0 x 10 4 , 0.6 x 10 6 ), all in meters. 

Given Eq. (13) derived in Problem 3.3 and the values 
of the parameters A, B, and C, estimate the accelera- 
tion error due to the effect of the solar corona on the 
propagation of radio waves between the Earth and the 
spacecraft. 

Solution 3.4. The correction to the Doppler fre- 
quency shift is obtained from Eq. (13) by simple time 
differentiation. (The impact parameter depends on time 
as p = p(t) and may be expressed in terms of the rela- 
tive velocity of the spacecraft with respect to the Earth, 
v s=s 30km/s). Use Eq. (13) for the steady-state solar 
corona effect on the radio-wave propagation through the 



solar plasma. The effect of the solar corona is expected 
to be small on the Doppler frequency shift - our main 
observable. The time-averaged effect of the corona on 
the propagation of the Pioneers' radio-signals is found to 
be small, of order 



±0.02 x 10~ 8 cm/s . 



(15) 



This small result is expected from the fact that most of 
the data used for the Pioneer analysis were taken with 
large Sun-Earth-spacecraft angles. 

Problem 3.5. The possibility that the Pioneer space- 
craft can hold a charge and be deflected in its trajec- 
tory by Lorentz forces was a concern due to the large 
magnetic field strengths near Jupiter and Saturn (sec 
Figs. 2 and l). 25 The magnetic field strength in the outer 
solar system is < 10~ 5 Gauss, which is about a factor 
of 10 5 times smaller than the magnetic field strengths 
measured by the Pioneers at their nearest approaches 
to Jupiter: 0.185 Gauss for Pioneer 10 and 1.135 Gauss 
for Pioneer 11. Furthermore, data from the Pioneer 10 
plasma analyzer can be interpreted as placing an upper 
bound of 0.1/uC on the positive charge during the Jupiter 
encounter. 26 

Estimate the upper limit of the contribution of the elec- 
tromagnetic force on the motion of the Pioneer spacecraft 
in the outer solar system. 

Solution 3.5. We use the standard relation, F = qv x 
B, and find that the greatest force would be on Pioneer 11 
during its closest approach, < 20 x 10~ 8 cm/s 2 . However, 
once the spacecraft reached the interplanetary medium, 
this force decreased to 



CTLorontz < 2 X 10 12 CHl/s 2 



(16) 



which is negligible. 

Problem 3.6. Long-term frequency stability tests are 
conducted with the exciter /transmitter subsystems and 
the Deep Space Network's (DSN) 27 radio-science open- 
loop subsystem. An uplink signal generated by the ex- 
citer is converted at the spacecraft antenna to a downlink 
frequency. The downlink signal is then passed through 
the RF-IF down-converter present at the DSN antenna 
and into the radio science receiver chain. This tech- 
nique allows the processes to be synchronized in the DSN 
complex based on the frequency standards whose Allan 
deviations are the order of a y <~ 10~ 14 -10 -15 for inte- 
gration times in the range of 10 s to 10 3 s. (The Allan 
deviation is the variance in the relative frequency, de- 
fined as y = Av/v, or a y ~ a v jv. See the discussion 
in Appendix A.) For the S-band frequencies of the Pi- 
oneers, the corresponding variances are 1.3 x 10~ 12 and 
1.0 x 10~ 12 , respectively, for a 10 3 s Doppler integration 
time. 

The influence of the clock stability on the detected 
acceleration, ap, can be estimated from the reported 
variances for the clocks, a y . The standard single mea- 
surement error on acceleration is derived from the time 
derivative of the Doppler frequency data as (c<j v )/t, 
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where the variance, a y , is calculated for 1000 s Doppler 
integration time and r is the signal averaging time. This 
relation provides a good rule of thumb when the Doppler 
power spectral density function obeys a 1// flicker noise 
law, which is approximately the case when plasma noise 
dominates the overall error in the Doppler observable. 

Assume a worst case scenario, where only one clock was 
used for the entire 11.5 years of data. Estimate the in- 
fluence of that one clock on the reported accuracy of the 
detected anomaly, ap. Combining a y — Au/v , the frac- 
tional Doppler frequency shift for a reference frequency of 
vq rs 2.295 GHz, with the estimate for the variance, a y = 
1.3 x 10 -12 , yields the upper limit for the frequency uncer- 
tainty introduced in a single measurement by the insta- 
bilities in the atomic clock: a v = vn<J y = 2.98 x 10~ 3 Hz 
for a 10 3 Doppler integration time. 

Obtain an estimate for the total effect of phase and fre- 
quency stability of clocks to the Pioneer anomaly value. 
How significant is this effect? 

Solution 3.6. To obtain an estimate for the total 
effect, recall that the Doppler observation technique is 
essentially a continuous count of the total number of 
complete frequency cycles during the observational time. 
The Pioneer measurements were made with duration 
w 10 3 s. Therefore, within a year we could have as many 
as N w 3, 150 independent single measurements of the 
clock. The values for a y and N yield an upper limit for 
the contribution of the atomic clock instability on the fre- 
quency drift: er c i oc k = o v j\fN w 5.3 x 10 _5 Hz/yr. The 
observed ap corresponds to a frequency drift of about 
0.2Hz/yr, which means that the error in ap is about 

(T freq = 0.0003 x 10~ 8 cm/s 2 , (17) 
which is negligible. 

IV. ON-BOARD SOURCES OF 
ACCELERATION NOISE 

In this section forces generated by on-board spacecraft 
systems that could contribute to the constant acceler- 
ation, ap, are considered. The on-board mechanisms 
discussed are the radio beam reaction force, the RTG 
heat reflecting off the spacecraft, the differential cmissiv- 
ity of the RTGs, the expelled Helium produced within 
the RTG, the thruster gas leakage, and the difference in 
experimental results from the two spacecraft. 

Problem 4.1. The Pioneers have a total nominal 
emitted radio power of 8W. The power is parametrized 
as 

P rp = / d6 smQV{6), (18) 
Jo 

where V(9) is the antenna power distribution. The radi- 
ated power is kept constant in time, independent of the 
coverage from ground stations. That is, the radio trans- 
mitter is always on, even when signals are not received 
by a ground station. 



The recoil from this emitted radiation produces an ac- 
celeration bias, 29 6 rp , on the spacecraft away from the 
Earth of 



where m is taken to be the Pioneer mass when half the 
fuel is gone, 241 kg, and j3 is the fractional component of 
the radiation momentum that is in the direction opposite 
to ap: 

1 ^max 

0=— d9 sin 9 cos 9V (9). (20) 

"rp Jo 

The properties of the Pioneer downlink antenna pat- 
terns are well known. The gain is given as (33.3 ±0.4) dB 
at zero (peak) degrees. The intensity is down by a factor 
of 2 (-3 dB) at 1.8°. It is down a factor of 10 (-10 dB) at 
2.7° and down by a factor of 100 (-20 dB) at 3.75°. (The 
first diffraction minimum is at a little over 4°.) There- 
fore, the pattern is a very good conical beam. 

Estimate the effect of the recoil force due to the emit- 
ted radio-power on the craft. What is the uncertainty in 
this estimation? 

Solution 4.1. Because cos 3.75° = 0.9978, we can take 
(3 = (0.99 ± 0.01), which yields 6 rp = 1.10. To estimate 
the uncertainty, we take the error for the nominal 8W 
power to be given by the 0.4 dB antenna error (0.10), 
and obtain 

a rp = 6 rp ± cr rp = (1.10 ± 0.10) x 10~ 8 cm/s 2 . (21) 




FIG. 3: A design picture of the SNAP 19 RTGs. Note that 
for better heat radiation, this final design had much larger 
fins than those in the earlier spacecraft concept of Fig. 2. 
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TABLE IV: Thermal system and on-board power for Pioneers 10 and 11. The core of the system is four SNAP-19 RTGs (SNAP 
stands for Space Nuclear Auxiliary Power). These RTGs are fueled with Plutonium 238 that undergoes nuclear a-decay via 
94PU 238 — ► 92 U 234 + 2He 4 , with a half life of 87.74 yr. The efficiency of the relevant heat-to-electric power conversion is about 
5 to 6%. 



Power available: 

> before launch electric total 165 W (by 2001 « 61 W) 

> needs 100 W to power all systems (e 24.3 W science instruments) 
Heat provided: 

D> before launch thermal fuel total 2580 W (by 2001 « 2050 W) 

> electric heaters; 12 one-W RHUs 

> heat from the instruments (dissipation of 70 to 120 W) 
Excess power/hcat: if electric power was > 100 W => 

t> thermally radiated into space by a shunt-resistor radiator, or 

> charge a battery in the equipment compartment 
Thermal control: to keep temperature within (0-90) F 

t> thermo-responsive louvers: completely shut/fully open for 40F/85F 
t> insulation: multi-layered aluminized mylar and kapton blankets 



Problem 4.2. It has been argued that the anomalous 
acceleration is due to anisotropic heat reflection off of the 
back of the spacecraft high-gain antennas, the heat com- 
ing from the RTGs. Before launch, the four RTGs had 
a total thermal fuel inventory of 2580 W (w 2050 W in 
2001). They produced a total electrical power of 165W 
(w 61 W in 2001). (In Table IV a description of the over- 
all power system is given.) Thus, by 2002 approximately 
2000 W of RTG heat still had to be dissipated. Because 
only ss 63 W of directed power could have explained the 
anomaly, in principle, there was enough power to explain 
the anomaly. 

The main bodies of the RTGs are cylinders and are 
grouped in two packages of two. Each package has two 
cylinders, end to end, extending away from the antenna. 
Each RTG has six fins separated by equal angles of 60° 
that go radially out from the cylinder. The approximate 
dimensions for the RTGs are: the cylinder and fin lengths 
are 28.2 and 26.0 cm, the radius of the central RTG cylin- 
der is 8.32 cm, and the RTG fin width as measured from 
the surface of the cylinder to the outer fin tip is 17.0 cm 
(see Fig. 3). 

The RTGs are located at the end of booms, and rotate 
about the spacecraft in a plane that contains the approx- 
imate base of the antenna. From the closest axial center 
point of the RTGs, the antenna is seen nearly "edge on" 
(the longitudinal angular width is 24.5°). The total solid 
angle subtended is ~ 1-2% of 4ir steradians. A more 
detailed calculation yields a value of 1.5%. 

Estimate the contribution of the RTG heat reflecting 
off the spacecraft to the Pioneer anomaly. Can you ex- 
plain the observed anomaly with this mechanism? Use 
the spacecraft geometry and the resultant RTG radiation 
pattern. 

Solution 4.2. There are two reasons that preclude 
this mechanism. One is the spacecraft geometry. Even if 
we take the higher bound of 2% of 4-7T steradians solid an- 
gle as the fraction of solid angle covered by the antenna, 



the equivalent fraction of the RTG heat could provide at 
most w 40 W. The second reason is the RTGs' radiation 
pattern. Our estimate was based on the assumption that 
the RTGs are spherical black bodies, which they are not. 

In fact, the fins are "edge on" to the antenna (the fins 
point perpendicular to the cylinder axes). The largest 
opening angle of the fins is seen only by the narrow angle 
parts of the antenna's outer edges. If we ignore these 
edge effects, only w 2.5% of the surface area of the RTGs 
is facing the antenna, which is a factor of 10 less than 
that obtained from integrating the directional intensity 
from a hemisphere: (f hcml d£l cos 9) / (4ir) = 1/4. Thus, 
we have only 4 W of directed power. 

The force from 4 W of directed power suggests a sys- 
tematic bias of w 0.55 x 10~ 8 cm/s 2 . If we add an un- 
certainty of the same size, we obtain a contribution from 
heat reflection of 

a hr = (-0.55 ± 0.55) x 10~ 8 cm/s 2 . (22) 

If this mechanism were the cause, ultimately an un- 
ambiguous decrease in the size of ap should be observed, 
because the RTGs' radioactively produced radiant heat 
is decreasing. As we noted, the heat produced is now 
about 80% of the original magnitude. In fact, we would 
expect a decrease of about 0.75 x 10~ 8 cm/s 2 in ap over 
the 11.5 year Pioneer 10 data interval if this mechanism 
were the origin of ap. 

Problem 4.3. Another suggestion related to the 
RTGs is the following: During the early parts of the 
missions, there might have been a differential change of 
the radiant emissivity of the solar-pointing sides of the 
RTGs with respect to the deep-space facing sides. Note 
that, especially closer in the Sun, the inner sides were 
subjected to the solar wind. On the other hand, the 
outer sides were sweeping through the solar system dust 
cloud. These two processes could result in different lev- 
els of degradation of the optical surfaces of the RTGs. In 
turn this degradation could result in asymmetric patterns 
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of heat radiation away from the RTGs in the more/aft di- 
rections along the spin axis. Therefore, it can be argued 
that such an anisotropy may have caused the anomaly. 

The six fins of each RTG, designed to provide the bulk 
of the heat rejection capacity are fabricated of mag- 
nesium alloy. 2 The metal is coated with two to three 
mils of zirconia in a sodium silicate binder to provide 
a high emissivity (w 0.9) and low absorptivity (rs 0.2) 
(see Fig. 3). 

Estimate the possible difference in the fore-and-aft 
emissivities of the RTGs needed to support this mech- 
anism. Can this mechanism explain the observed 
anomaly? Discuss the significance of radioactive decay 
for this mechanism. 

Solution 4.3. Depending on how symmetrically fore- 
and-aft they radiated, the relative fore-and-aft emissivity 
of the alloy would have had to have changed by w 10% to 
account for ap. Given our knowledge of the solar wind 
and the interplanetary dust, this amount of a radiant 
change would be difficult to explain, even if it had the 
right sign. 

To obtain a reasonable estimate of the uncertainty, 
consider if one side (fore or aft) of the RTGs has its 
emissivity changed by only 1% with respect to the other 
side. In a simple cylindrical model of the RTGs, with 
2000 W power (only radial emission is assumed with no 
loss out of the sides), the ratio of the power emitted by 
the two sides would be 995/1005 = 0.99, or a differen- 
tial emission between the half cylinders of 10 W. There- 
fore, the fore/aft asymmetry toward the normal would 
be 10W x i f* d<(> sin 0« 6.37 W. 

A more sophisticated model of the fin structure starts 
off from ground zero to determine if we still obtain about 
the order of magnitude. 2 Indeed we do, 6.12W, which is 
the value we use in the following. (We refer the reader 
to Ref. 2 for this discussion, because the basic physics 
already has been discussed.) 

If we take 6.12 W as the uncertainty from the differ- 
ential emissivity of the RTGs, we obtain an acceleration 
uncertainty of 

o-dc = 0.85 x lO^cm/s 2 . (23) 

Note that lOcdc almost equals our final result for ap 
(see Section VI). This correspondence is the origin of our 
previous statement that w 10% differential emissivity (in 
the correct direction) would be needed to explain ap. 

Finally, consider the significance of radioactive decay 
for this mechanism. The formal statistical error of the 
determination (before considering experimental system- 
atics) was small 2 (see Sec. V). An example is a specific 
one-day "batch-sequential filtering" 28 value for ap, av- 
eraged over the entire 11.5 year interval. This value was 
a P = (7.77±0.16)xl0~ 8 cm/s 2 , where a batch-sequential 
error tends to be large. From radioactive decay, the value 
of ap should have decreased by 0.75 x 10~ 8 cm/s 2 over 
11.5 years, which is five times the large batch sequential 
variance. Even more stringently, this bound is good for 
all radioactive heat sources. So, what if we were to argue 



that emissivity changes occurring before 1987 were the 
cause of the Pioneer effect? There still should have been 
a decrease in ap with time since then, which has not been 
observed. 

Problem 4.4. Another possible on-board systematic 
error is from the expulsion of the He being created in the 
RTGs from the a-decay of 238 Pu. 

The Pioneer RTGs were designed so that the He pres- 
sure is not totally contained within the Pioneer heat 
source over the life of the RTGs. Instead, the Pioneer 
heat source contains a pressure relief device that allows 
the generated He to vent out of the heat source and 
into the thermoelectric converter. The thermoelectric 
converter housing-to-power output receptacle interface is 
sealed with a viton O-ring. The O-ring allows the helium 
gas within the converter to be released by permeation 
into the space environment throughout the mission life 
of the Pioneer RTGs. 

The actual fuel (see Fig. 3) is formed into a stack of 
disks within a container, each disk having a height of 
0.212" and a diameter of 2.145". With 18 disks in each 
RTG and four RTGs per mission, the total volume of fuel 
is about 904 cm 3 . The fuel is a plutonium moly cermet 
and the amount of 238 Pu in this fuel is about 5.8 kg. With 
a half life of 87.74 yr, the rate of He production (from 
plutonium a-decay) is about 0.77 g/yr, assuming it all 
leaves the cermet and goes into the RTG chamber. To 
make this mechanism work, He leakage from the RTGs 
must be preferentially directed away from the Sun with 
a velocity large enough to cause the acceleration. 

Can this mechanism explain the Pioneer anomaly? 
What is the largest effect possible with this mechanism? 
What is the uncertainty? 

Solution 4.4. The operational temperature on the 
RTG surface at 433 K implies a helium velocity of 
1.22 km/s. We use this value in the rocket equation, 

a(t) = -v(t)j t [lnM(t)}. (24) 

We assume that the Pioneer mass corresponds to half of 
the fuel and that the gas is all unidirected and find a 
maximum bound on the possible acceleration of 1.16 x 
10~ 8 cm/s 2 . Thus, we can rule out helium permeating 
through the O-rings as the cause of ap. 

If we assume a single elastic reflection, we can esti- 
mate the acceleration bias 29 in the direction away from 
the Sun. This estimate is (3/4) sin 30° times the average 
of the heat momentum component parallel to the short- 
est distance to the RTG fin. By using this geometric in- 
formation, we find the bias would be 0.31 x 10~ 8 cm/s 2 . 
This bias effectively increases the value of the solution for 
ap, which is too optimistic given all the complications of 
the real system. Therefore, we can take the systematic 
expulsion to be 

a Hc = (0.15 ± 0.16) x 10~ 8 cm/s 2 . (25) 

Problem 4.5. The effect of propulsive mass expul- 
sion due to gas leakage has to be assessed. Although 
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this effect is largely unpredictable, many spacecraft have 
experienced gas leaks producing accelerations on the or- 
der of 10~ 7 cm/s 2 . Gas leaks generally behave differently 
after each maneuver. The leakage often decreases with 
time and becomes negligible. 

Gas leaks can originate from Pioneer's propulsion sys- 
tem, which is used for mid-course trajectory maneu- 
vers, for spinning up or down the spacecraft, and for 
orientation of the spinning spacecraft. The Pioneers 
are equipped with three pairs of hydrazine thrusters, 
which are mounted on the circumference of the Earth- 
pointing high gain antenna. Each pair of thrusters forms 
a thruster cluster assembly with two nozzles aligned op- 
posite to each other. For attitude control, two pairs of 
thrusters can be fired forward or aft and are used to pre- 
cess the spinning antenna. The other pair of thrusters 
is aligned parallel to the rim of the antenna with noz- 
zles oriented in co- and contra-rotation directions for 
spin/despin maneuvers. 

During both observation intervals for the two Pio- 
neers, there were no trajectory or spin/despin maneu- 
vers. Thus, in this analysis we are mainly concerned 
with precession (that is, orientation or attitude control) 
maneuvers only. Because the valve seals in the thrusters 
can never be perfect, we ask if the leakages through the 
hydrazine thrusters could be the cause of the anomalous 
acceleration, ap. 

Consider the possible action of gas leaks originating 
from the spin/despin thruster cluster assembly. Each 
nozzle from this pair of thrusters is subject to a certain 
amount of gas leakage. But only a differential leakage 
from the two nozzles would produce an observable effect, 
causing the spacecraft to either spin down or up. 

To obtain a gas leak uncertainty (and we emphasize un- 
certainty rather than "error" because there is no other ev- 
idence), we ask, "How large a differential force is needed 
to cause the spin down or spin up effects that are ob- 
served?" The answer depends on the moment of inertia 
about the spin axis, I z s=s 588.3 kg m 2 , the antenna ra- 
dius, R = 1.37m, and the observed spin down rates over 
three time intervals. In Rcf. 2 the final effect was deter- 
mined at the end by considering the spin-downs over the 
three time periods to be manifestations of errors that are 
uncorrelated and are normally distributed around zero 
mean. For simplicity, we include here the effect of this 
analysis in terms of an overall effective spin-down rate of 
(9 = -0.0067 rpm/yr. 

Estimate the uncertainty in the value for the Pioneer 
anomaly from the possibility of undetected gas leaks. 

Solution 4.5. If we use the information in Problem 
4.5 and take the antenna radius, R, as the lever arm, we 
can calculate the differential force needed to torque the 
spin-rate change: 

T Q 

F § = ^- = 0.95 x 10~ 3 dyn. (26) 
R 

It is possible that a similar mechanism of undetected 
gas leakage could be responsible for the net differential 



force acting in the direction along the line of sight. In 
other words, what if there were some undetected gas leak- 
age from the thrusters oriented along the spin axis of the 
spacecraft that is causing a pi How large would this leak- 
age have to be? A force of 

F ap = map = 21.11 x 10~ 3 dyn (27) 

(m = 241 kg) would be needed to produce our final unbi- 
ased value of ap. Given the small amount of information, 
we conservatively take as our gas leak uncertainties the 
acceleration values that would be produced by differen- 
tial forces equal to 

F g] w ±V2F § = ±1.35 x 10~ 3 dyn. (28) 

The argument for the result obtained in Eq. (28) is that 
we are accounting for the differential leakages from the 
two pairs of thrusters with their nozzles oriented along 
the line of sight direction. Equation (28) directly trans- 
lates into the acceleration errors introduced by the leak- 
age during the interval of Pioneer 10 data, 

cr g i = ±F g i/m = ±0.56 x 10~ 8 cm/s 2 . (29) 

At this point, we must conclude that the gas leak mech- 
anism for explaining the anomalous acceleration is very 
unlikely, because it is difficult to understand why it would 
affect Pioneers 10 and 11 by the same amount. We also 
expect a gas leak would contribute to the anomalous ac- 
celeration, ap, as a stochastic variable obeying a Poisson 
distribution. Instead, analyses of different data sets in- 
dicate that a p behaves as a constant bias rather than as 
a random variable. 2 

Problem 4.6. There are two experimental results for 
the Pioneer anomaly from the two spacecraft: 7.84 x 
10~ 8 cm/s 2 (Pioneer 10) and 8.55 x 10~ 8 cm/s 2 (Pioneer 
11). The first result represents the entire 11.5 year data 
period for Pioneer 10; Pioneer ll's result represents a 
3.75 year data period. 

The difference between the two craft could be due to 
different gas leakage. But it also could be due to heat 
emitted from the RTGs. In particular, the two sets of 
RTGs have had different histories and so might have dif- 
ferent emissivities. Pioneer 11 spent more time in the 
inner solar system (absorbing radiation). Pioneer 10 has 
swept out more dust in deep space. Further, Pioneer 11 
experienced about twice as much Jupiter/Saturn radia- 
tion as Pioneer 10. However, if the Pioneer effect is real, 
and not due to an extraneous systematic like heat, these 
numbers should be approximately equal. 

Estimate the value for the Pioneer anomaly based on 
these two independent determinations. What is the un- 
certainty in this estimation? 

Solution 4.6. We can calculate the time-weighted 
average of the experimental results from the two craft: 
[(11.5)(7.84) + (3.75)(8.55)]/(15.25) = 8.01 in units of 
10~ 8 cm/s 2 . This result implies a bias of o 2 , craft = 
0.17 x 10~ 8 cm/s 2 with respect to the Pioneer 10 exper- 
imental result ap( eX p) (see Eq. (35) below). In addition, 
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we can take this number to be a measure of the uncer- 
tainty from the separate spacecraft measurements, so the 
overall quantitative measure is 

0-2, craft = &2, craft ± 02, craft 

= (0.17 ±0.17) x l(T 8 cm/s 2 . (30) 

V. SOURCES OF COMPUTATIONAL ERRORS 

Given the very large number of observations for the 
same spacecraft, the error contribution from observa- 
tional noise is very small and not a meaningful measure of 
uncertainty. It is therefore necessary to consider several 
other effects in order to assign realistic errors. The first 
consideration is the statistical and numerical stability of 
the calculations. Then there is the cumulative influence 
of all program modeling errors and data editing deci- 
sions. Besides these factors, there are errors that may be 
attributed to the specific hardware used to run the orbit 
determination computer codes, together with the algo- 
rithms and statistical methods used to obtain the solu- 
tion. These factors have been included in the summary 
of the biases and uncertainties given in Table V. The 
largest contributor in this category, and the question we 
will discuss here, is the reason for and significance of a 
periodic term that appears in the data. 

Problem 5.1. In addition to the constant anomalous 
acceleration term, an annual sinusoid has been reported. 2 
The peaks of the sinusoid occur when the Earth is exactly 
on the opposite side of the Sun from the craft, where 
the Doppler noise is at a maximum. The amplitude of 
this oscillatory term by the end of the data interval was 
w 0.22 x 10~ 8 cm/s 2 . The integral of a sine wave in the 
acceleration, with angular velocity lu and amplitude do, 
yields the following first-order Doppler amplitude, in the 
fractional frequency change accumulated over the round- 
trip travel time of the radio-signal: 

^ = ^. (31) 

V CUJ 

The resulting Doppler amplitude from the sine wave with 
amplitude of cto = 0.22 x 10~ 8 cm/s 2 and annual angu- 
lar velocity u> = 2 x 10~ 7 rad/s is Av/is = 7.3 x 10~ 13 . 
At the Pioneer downlink S-band carrier frequency of 
w 2.29 GHz, the corresponding Doppler amplitude is 
0.002 Hz, that is, 0.11 mm/s. 

A four-parameter, nonlinear, weighted, least-squares 
fit to the annual sine wave was found with the amplitude 
v at = (0.1053 ±0.0107) mm/s, phase = (-5.3° ± 7.2°), 
angular velocity a> at = (0.0177 ± 0.0001) rad/day, and 
bias = (0.0720 ± 0.0082) mm/s. Standard data weight- 
ing procedures 2 yield post-fit weighted rms residuals of 
<jt = 0.1 mm/s (the Doppler error averaged over the data 
interval T). 

The amplitude, u at , and angular velocity w at , of the 
annual term result in a small acceleration amplitude of 
flat = Wat^at = (0.215 ± 0.022) x lO^cm/s 2 . The cause 



is most likely due to errors in the navigation programs' 
determinations of the direction of the spacecraft's orbital 
inclination to the ecliptic. Estimate the annual contribu- 
tion to the error budget for ap. 

Solution 5.1. First observe that the standard errors 
for the radial velocity v r , and acceleration, a r , are es- 
sentially what would be expected for a linear regression. 
The caveat is that they are scaled by the root sum of 
squares of the Doppler error and unmodeled sinusoidal 
errors, rather than just the Doppler error. Furthermore, 
because the error is systematic, it is unrealistic to assume 
that the errors for v r and a r can be reduced by a factor 
1/y/N, where N is the number of data points. Instead, if 
we average their correlation matrix over the data interval, 
T, we obtain the estimated systematic error of 

< = ^< = ^(4 + < t ), (32) 

where ot — 0.1 mm/s is the Doppler error averaged over 
T (not the standard error of a single Doppler measure- 
ment), and <7„ at is equal to the amplitude of the unmod- 
eled annual sine wave divided by The resulting root 
sum of squares error in the radial velocity determination 
is about a Vr = (erf, + cr 2 ) 1 / 2 = 0.15mm/s for both Pio- 
neer 10 and 11. Our values of ap were determined over 
time intervals of longer than a year. At the same time, 
to detect an annual signature in the residuals, we need at 
least half of the Earth's orbit to be complete. Therefore, 
with T = 1/2 yr, Eq. (32) results in an acceleration error 
of 

0.50 mm/s „ « ,2 ,„„x 

cr ar = ^—!— = 0.32 x 10~ 8 cm/s 2 . (33) 

This number is assumed to be the systematic error from 
the annual term. 



VI. SUMMARY OF ERRORS AND FINAL 
RESULT 

The tests discussed in the preceding sections have 
considered various potential sources of systematic er- 
ror. The results are summarized in Table V. The only 
major source of uncertainty not discussed here is the 
non-isotropic radiative cooling of the spacecraft. 2,9 Other 
more minor effects are discussed in Ref. 2. Table V 
summarizes the systematic errors affecting the measured 
anomalous signal, the "error budget" of the anomaly. It 
is useful for evaluating the accuracy of the solution for ap 
and for guiding possible future efforts with other space- 
craft. 

Finally, there is the intractable mathematical problem 
of how to handle combined experimental systematic and 
computational errors. In the end it was decided to treat 
them all in a least squares uncorrelated manner. 2 

Problem 6.1. The first column in Table V gives the 
bias, &p, and the second gives the uncertainty ±op. The 
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TABLE V: Error budget: A summary of biases and uncertainties. 



Item Description of error budget constituents 


Bias 

10 -8 cm/s 2 


Uncertainty 
10~ 8 cm/s 2 


1. 


Sources of acceleration noise external to the craft 








(a) Solar radiation pressure 


+0.03 


±0.01 




(b) Solar wind 




± < 10~ 3 




(c) Solar corona 




±0.02 




(d) Electromagnetic Lorentz forces 




± < io- 4 




(e) Phase stability and clocks 




± < 0.001 


2. 


On-board sources of acceleration noise 








(a) Radio beam reaction force 


+ 1.10 


±0.10 




(b) RTG heat reflected off the craft 


-0.55 


±0.55 




(c) Differential emissivity of the RTGs 




±0.85 




(d) Non-isotropic radiative cooling of the spacecraft 




±0.48 




(e) Expelled Helium produced within the RTGs 


+0.15 


±0.16 




(f ) Gas leakage 




±U.OD 




(g) Variation between spacecraft determinations 


+0.17 


±0.17 


3. 


Sources of computational errors 








(a) Accuracy of consistency /model tests 




±0.13 




(b) Annual term 




±0.32 




(c) Other important effects 




±0.04 




Estimate of total bias/error 


+0.90 


±1.33 



RTGs, radiative cooling, and gas leaks, ±0.85, ±0.48, 
and ±0.56, respectively, xl0~ 8 cm/s 2 . 

The computational errors are listed in the third group 
of Table V. We obtain b P ± a P = 0.90 x 10~ 8 cm/s 2 and 
ap = 1.33 x 10~ 8 cm/s 2 . Therefore, from Eq. (34) the 
final value for ap is 

a P = (8.74 ± 1.33) x 10~ 8 cm/s 2 . (36) 

The effect is clearly significant and remains to be ex- 
plained. 
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constituents of the error budget are listed in three cate- 
gories: systematic errors generated external to the space- 
craft; on-board generated systematic errors, and compu- 
tational errors. Our final result will be an average 

ap = ap(o X p) +bp±cr P , (34) 

where 

ap(cxp) = (7.84 ± 0.01) x 10" 8 cm/s 2 (35) 

is our formal solution for the Pioneer anomaly that was 
obtained with the available data set. 2 

Discuss the contributions of various effects to the total 
error budget and determine the final value for the Pioneer 
anomaly, ap. 

Solution 6.1. The least significant factors of our error 
budget are in the first group of effects, those external to 
the spacecraft. From Table V we see that some are near 
the limit of contributing, and in total they are insignifi- 
cant. 

The on-board generated systematics are the largest 
contributors to the total. All the important constituents 
are listed in the second group of effects in Table V. 
Among these effects, the radio beam reaction force pro- 
duces the largest bias, 1.10 x 10 _8 cm/s 2 , and makes 
the Pioneer effect larger. The largest bias/uncertainty 
is from RTG heat reflecting off the spacecraft. The effect 
is as large as (—0.55 ± 0.55) x 10~ 8 cm/s 2 . Large un- 
certainties also come from differential emissivity of the 
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APPENDIX A: ALLAN DEVIATION deviation, which defined as 



To achieve the highest accuracy, the Doppler frequency 
shift of the carrier wave is referenced to ground-based 
atomic frequency standards. The observable is the frac- 
tional frequency shift of the stable and coherent two-way 
radio signal (Earth-spacecraft-Earth) 

where and vt are, respectively, the transmitted and 
received frequencies, t is the receiving time, and 2L is 
the overall optical distance (including diffraction effects) 
traversed by the photon in both directions. 

The received signal is compared with the expected 
measurement noise in y, given by Eq. (Al). The qual- 
ity of such a measurement is characterized by its Allan 



^(r) = ^<(w+i(r)-y i (r))2) ) (A2) 

where yi — \ J^ +T y(t)dt. The Allan deviation is the 
most widely used figure of merit for the characterization 
of frequency in this context. 

For advanced planetary missions such as Cassini, 
equipped with multi-frequency links in the X- and Ka- 
bands, Allan deviations reach the level of a y « 10~ 14 
for averaging times between 10 3 and 10 4 s. The Pioneer 
spacecraft were equipped with S-band Doppler commu- 
nication systems for which Allan deviations were usually 
the order ofl of 1.3 x 10~ 12 and 1.0 x 10~ 12 , respectively, 
for 10 3 s Doppler integration times. 
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To determine the state of the spacecraft at any given time, 
the Jet Propulsion Laboratory's orbit determination pro- 
gram uses a batch-sequential filtering and smoothing algo- 
rithm with process noise. 2 For this algorithm, any small 
anomalous forces may be treated as stochastic parameters 
affecting the spacecraft trajectory. As such, these parame- 
ters are also responsible for the stochastic noise in the ob- 
servational data. To better characterize these noise sources, 
the data interval is split into a number of constant or 
variable size intervals, called "batches," and assumptions 
are made for statistical properties of these noise factors. 
The mean values and the second moments of the unknown 
parameters are estimated within the batch. In batch se- 
quential filtering, the reference trajectory is updated after 
each batch to reflect the best estimate of the true trajec- 
tory. (For more details consult Ref. 2 and also references 
therein.) 

We use the symbol b to identify biases in the measured 
value of the Pioneer anomaly produced by various sources 
of acceleration noise. The uncertainly in the contribution 
of these noise sources to the anomalous acceleration is con- 
ventionally denoted by a. See an example of this usage in 
Eq. (21). 



